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Abstract. We consider the Dirichlet Laplacian in a waveguide of uniform width and infinite 
length which is ideally divided into three parts: a "vertex region", compactly supported and 
■ with non zero curvature, and two "edge regions" which are semi-infinite straight strips. We 

make the waveguide collapse onto a graph by squeezing the edge regions to half-lines and the 
pL^ , vertex region to a point. In a setting in which the ratio between the width of the waveguide 

1 and the longitudinal extension of the vertex region goes to zero, we prove the convergence of 

' the operator to a selfadjoint realization of the Laplacian on a two edged graph. In the limit 

operator, the boundary conditions in the vertex depend on the spectral properties of an effective 
i i ' one dimensional Hamiltonian associated to the vertex region. 

O |' Keywords: Quantum graphs, quantum waveguides. 

! MSG 2010: 81Q35, 81Q37, 47A10. 

1. Introduction 

i— i' Metric graphs are objects of great interest as simple but highly not trivial tools to investigate 

^ . questions covering several topical issues in pure and applied mathematics. They offer a twofold 

advantage: they are quite simple objects (in many cases formulas for the relevant quantities can 
C^~- ■ be explicitly computed), showing at the same time non trivial features. We evade the duty of 

giving an exhaustive list of applications of metric graphs and limit ourselves to point out two 
recent volumes on the subject [31 where an extensive list of references can be found. 
Our work is addressed to investigate a long standing problem within the theory of metric graphs, 
which is finding a rigorous justification for the use of graph models to approximate dynamics 
in networks of thin waveguides. In many applications the problem is reduced to compare, 
in some suitable sense, the Laplacian in a squeezing network with a selfadjoint realization of 
£> | the Laplacian on the graph. A difficult feature of the problem originates from the fact that 

requiring selfadjointness does not fix univocally the definition of the operator on the graph. 
Each selfadjoint realization of the Laplacian on a graph is associated to certain boundary (or 
matching) conditions in the vertices. There are different ways to express the most general 
selfadjoint boundary conditions, see, e.g., [8j and [25], for the sake of simplicity let us specify 
one of them for the case of a star- graph Q, that is a graph with N edges of infinite length and 
one vertex v. We denote by Mg the Hilbert space naturally associated to Q, a function x £ Mg 
is a vector with N components, x = (x±, ...,xn), such that Xj S L 2 ((0, oo)). In each edge the 
vertex v is identified with the origin. Consider an orthogonal projector n acting in C^, and 
a selfadjoint operator O in Ran(n _L ) C C^, with H 1 - = 1 — 11. We denote by — Ag' the 
selfadjoint realization of the Laplacian in Mg associated to the couple of operators (n, B). Let 
x = (xi, ...,xn) be a vector in the domain of —Ag'~, then Xj £ H 2 ((0,oo)) and x must satisfy 
the boundary conditions 

Ux(v) = , uV(u) + QU^xiv) = , 
where we denoted by x(v) and x'(v) the vectors in defined by x(v) = (xi(0), iejv(O)) and 



X 



x'(v) = (2^(0), x^(0)) , moreover — Ag ,e x = ( 



1 



2 



Among the possible realizations of —Ag'~ we mention: 

- The Dirichlet (or decoupling) Laplacian: it is defined by the boundary condition x(v) = 
0, which corresponds to the choice II = 1, in this case dim[Ran (IT 1 -)] = 0. 

- The weighted Kirchhoff Laplacian: it is defined by the conditions ajXi(0) = ajXj(O), for 
all i 7^ j, and c^ x 'i(0) = 0, with cti G C. This conditions correspond to the choice 

(1-1) (1% = Sij - - , i,j = l,...,N, 8 = 0, 

where 5{j denotes the Kronecker symbol. In this case dim[Ran (U"- 1 -)] = 1. When all the 
constants ctj are equal the operator defined by the choice (jl.ip is the so called Kirchhoff 
(or standard) Laplacian and the boundary conditions read £i(0) = ... = Xjv(0) and 

EiI a ^(o) = o. 

A central problem in using operators in the family — A^' to approximate the dynamics gener- 
ated by the Laplacian in a network of thin tubes squeezing to a graph is to understand which 
boundary conditions in the vertex arise in limit (see, e.g., [13] for a review on this topic). It 
comes out that this problem strongly depends on what kind of Laplacian is taken in the squeezing 
network. 

With Neumann conditions on the boundary of the squeezing network (or if one considers a 
network without boundary see, e.g., [141 130]) the problem is in some sense easier. Under general 
assumptions on the network properties (see the discussion in [2]), the limit operator on the 
graph is characterized by weighted Kirchhoff conditions in the vertex, with ctj real and the ratio 
oti/ctj being related to the relative radius of network's edges. 

At the best of our knowledge the first results in a setting of Neumann networks concerned 
the analysis of the convergence of stochastics processes associated with elliptic operators on 
the network, see [18] (see also [2], for recent results on the analysis of stochastics processes in 
networks of thin tubes). In [26, 33KH] the convergence of eigenvalues was analyzed while in [30] 
it was proved the convergence of operators in norm resolvent sense (see also [Ml EH] ) . We also 
mention the work [5] in which the problem is analyzed by using the theory of Gromov-Hausdorff 
convergence of manifolds and the techniques developed in [27] , 

In the setting with Neumann boundary conditions we quote also: the work [15] in which the 
convergence of resonances is analyzed and the work [16] where the presence of a scaled potential 
supported in the vertex region is taken into account. It is also worth pointing out that there 
exist works addressed to the analysis of the corresponding nonlinear problem, see, e.g., [31] , and 
references therein, for the case of squeezing tubular domains and [24J in which it is studied the 
convergence of semilinear elliptic equations. 

The case of a Dirichlet network squeezing to a graph is more complicated and not yet completely 
understood. One main difficulty is related to the fact that as the width of the network squeezes 
to zero the energy functional diverges as the inverse of the square of the width. Then to get a 
meaningful limit a renormalization of the energy is needed. Moreover geometrical perturbations 
can modify substantially the spectral properties of the Laplacian in the network, a well known 
phenomenon is, for example, the appearing of isolated eigenvalues in curved waveguides, see 
[17] . Similar problems arise when the confinement is obtained by an holonomic constrain, see, 

e.g., m- 

In |29] it is shown that in a Dirichlet network, if the volume of the vertex region is small enough, 
then the boundary conditions in the vertex are of Dirichlet type. Nevertheless it is known 
that more general assumptions cannot exclude different limits. In particular it is understood 
that weighted Kirchhoff conditions of the form described above may arise whenever an effec- 
tive Hamiltonian associated to the vertex region exhibits eigenvalues at the threshold of the 
transverse energy see, e.g., |20 t 1211 [28] and the recent work [10] in which these ideas are further 
exploited and the occurrence of non homogenous terms in the weighted Kirchhoff condition is 
pointed out. 
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In a previous work pQ we considered the Laplacian with Dirichlet conditions on the boundary 
of a planar waveguide of uniform width collapsing onto a two edged graph. We proved that two 
possible limits can arise: 

- In one case, which we call generic, the limit operator on the graph is defined by the 
Dirichlet condition in the vertex. 

- In a second case, which we call non-generic, the limit operator on the graph is defined 
by a weighted Kirchhoff condition in the vertex. The occurrence of the non-generic case 
is related to the existence of a zero energy resonance for an effective one dimensional 
Hamiltonian associated with the squeezing waveguide. 

A crucial hypothesis within the model analyzed in [Ij is the assumption that exist two scales of 
energy: a short one associated to the width of the waveguide, 5; a large one associated to the 
longitudinal extension of the vertex region, e, with 5 < e 5//2 . This assumption implies that the 
dynamics in the waveguide is adiabatically separated into a "fast" component associated to the 
transverse coordinate and a "slow" component associated to the longitudinal coordinate, making 
possible to analyze the problem in two steps. First the dynamics is reduced to a one dimensional 
problem by a projection onto the transverse eigenfunctions and taking the limit 5 — > (in norm 
resolvent sense); then the operator on the graph can be obtained by taking the limit e — > (in 
norm resolvent sense) of the effective one dimensional Hamiltonian. For the analysis of a similar 
problem in three dimensional waveguides with methods of Gamma-convergence we mention the 
paper jS], see also [S]. 

The weak point of this model, and of the approach used in [T], is that it does not offer enough 
margin for generalizations to settings in which the underlying graph has three or more edges. 
Aim of the present work is to give a better understanding of the result obtained in [I] in view of 
applications to multi-edged graphs. We review the problem separating explicitly the edges and 
vertex regions in the waveguide. According to the ideas in [TOj ES Ell EE], we reinterpret the 
existence of the zero energy resonance (which leads to coupling conditions in the vertex) as an 
eigenvalue at the energy of the transverse modes for an effective Hamiltonian associated to the 
vertex region. Moreover in the generic case we obtain better estimates (as compared to [I]) for 
the rate of convergence (5 < e 3//2 , instead of S < e 5 ^ 2 ). 

Related to the topic of our paper are also several works which concern the approximation of 
generic matching conditions in the vertex by scaled Schrodinger operators on the graph, see for 
example [8] and the recent review [13] . We also point out the paper [19] (and references therein) 
concerning the analysis of the convergence of Schrodinger operators in dimension one to <5'-type 
point interactions. 

Here is a short summary of our paper. In section [2] we define the waveguide as the union of 
three parts: the edges and vertex regions. Moreover we define the Dirichlet Laplacian in the 
waveguide. In section [3] we present our main results (Theorem [1] and Theorem [2]). Sections H] 
and [5] are devoted to the proofs of the main theorems. In section [6] we explain how to interpret 
the result stated in Theorems [1] and El in terms of the convergence (in norm resolvent sense) 
to an operator on the graph, see Theorem [3j We close the paper with section [7] in which we 
summarize our results and point out some generalizations. 



2. The model 

We consider a smooth waveguide of uniform width 5 > embedded in R 2 . We assume that the 
waveguide can be decomposed into three parts: two straight edges Ei t s and I?2,<5 and one vertex 
region V$ )E , e > being a parameter characterizing the "longitudinal extension" of the vertex 
region. Throughout the paper we shall always assume that 5/e ^ 1. Each edge Ej ; g, j = 1,2, 
can be identified with the manifold (E, f)^) where E := (0, oo) x (0, 1) and t)$ is the metric 



fa ■= ds 2 + 5 2 du 2 , 



where s G (0, oo) and u G (0, 1). The vertex region can be identified with the manifold (V,Qs,e) 
where V := (— 1, 1) x (0, 1) and Qs t£ is the metric 

qs, £ ■■=e 2 g S / £ ds 2 + 5 2 du 2 

with s G (—1,1) and u G (0, 1). The function g$/ £ (s,u) is defined by 

(2.1) g s/£ (s,u) = (l + u5/e 1 (s)) 2 , 

where 7(s) is a function of s and we assume that 7 G C^°((— 1,1)) and ||'y[|i 00 ((— 1,1)) < 1- We 
notice that this implies that gs/ £ G C 0O (l A ) and that for all < S ^ e ^ 1, the bounds 

( 2 - 2 ) < (1 - ||7||l-»((-i,i))) 2 < \\95/s\\l<x>{V) ^ (l + I|7||l-((-i,i))) 2 < 4 

hold true. The waveguide is obtained by identifying the boundary of E\j corresponding to s = 

with the boundary of Vs ;£ corresponding to s = — 1 and the boundary of E 2 s corresponding to 

s = with the boundary of Vs :£ corresponding to s = 1. 

We denote by Ms :£ the complex Hilbert space 

M SjE := L 2 (£ li5 ) L 2 (E 2 , S ) L 2 (^, e ) , 

where L 2 (E j: s) = L 2 (-E,f)a), j = 1,2, and L 2 (Va i£ ) = L 2 (V, gs,e)- Moreover we denote by a 
capital Greek letter a vector in Ms ;£ , 9 = (V>i, V>2> VvO- Given two vectors 3> = (<j>i, (j) 2 , <j> v ) G Ms 
and ^ = (ipi,ip2,ipv) G Ha )e the scalar product (<£, ^)H a g and the norm H^Hhj e read respectively 

roo rl rl rl 

($,V) Mse = V / / ^Vfcdet[^] 1/2 ds(iu+ / / ^^det[0,5 ie ] 1/2 dsdu, 

fc=12^° ^° J -I JO 

ll*lk, e = [(*>* 

with ~~ denoting the complex conjugation and 
(2.3) det[fj a ] = 5 2 , det[fla,J = 5 2 e 2 5 5 / £ 

In Ma e we define the sesquilinear form Qs :£ 

f 00 rl 



£ 



1 1/2 



n , T ,i /"°° / \ d( t>kd^k . 1 dfadipk] 1/2 



+ 



1 rl 



1 JO 



1 d(f) v dip v 1 <9</>„ dip v 



detlQse^^dsdu . 



e2 98/e ds ds 5 2 du du 

Below we shall replace the term sesquilinear form by quadratic form (this is justified by the 
polarization theorem, see, e.g., [23]) and we shall also use the notation (5a, £ [ v I'] = Qs )£ [^f, Let 
C°° be the set 
(2.4) 

C°° :={* = (^1,^2,^)1^1,^2 G C^°(E), V, GC°°(y); 

^1(5,0) = ^i(s,l) = V2M) = MsA) =Ms,0) = M*A) = 0; 

[d k s ^](0,u) = [(-e)- k dM(-hu); [d^ 2 ](0,u) = [ E - fc ^](l,it), Vk G N } , 

where we denoted by E the closure of E. The domain D(Qs t£ ) of the quadratic form Qs :£ is the 
closure of C°° equipped with the norm 

We denote by Hs j£ the unique selfadjoint operator in Ms j£ associated to the quadratic form Qs t£ 

(2.5) D(H S , £ ) := {* G D(Q 5 , £ )\ V$ G D(Q 5 , £ ) , Q a , e [$, *] = ($,H) H ,, e ; 3 G M 5 , £ } 

(2.6) H s ^ := 5 

(see, e.g., [23J for the concept of selfadjoint operators associated to closed quadratic forms). 



2.1. Unitarily equivalent Hamiltonian. Let L 2 (E) and L 2 (V) be the complex Hilbert spaces 
endowed with the norms 



OO fl 



U 3\\L 2 {E) 




\if}j\ 2 dsdu , j = 1,2; 



/o J o 

Let us denote by fi e the complex Hilbert space 

~2( Z7A ^ t2i Z7A rrs T 2 



i ./) 




Given two vectors <E> = 
and the norm II^IL read 



M £ := L (E) © L (E) © L (V, edsdu) . 
\.02; 0v) G fie an d ^ = (~ipi,i()2,ipv) G fie the scalar product ($,^) I fj 



(*.*)i.= E 



oo /■! . 



k=l,2 



1 /■! . 



4>kipkdsdu + e 



/0 ■/ 

|2 



1 JO 



cj) v ip v dsdu , 

v 1/2 



(2.7) 



II*IIh £ = (Jl^ilH^) + HtM^s) +e||^|l2 2( v)J 
For all < e ^ 1 we denote by C/<j )£ the unitary map Us <e '■ H^g — > H e defined by 

$i,Mv) = U S ^ U M V ) -(det^^Vi^etiy^^-^detfg^] 1 / 4 ^) 

=(5 1 /Vi,^ 1/2 ^,<5 1/2 ^) ) 
where we used equation (|2.3[) . 

For all < <5 ^ £ ^ 1, we denote by i?| £ the closure of C°°, defined in equation (|2.4p . with 
respect to the norm 
(2.8) 



11*11* : =n*nl s +E \\Mj\\h{Ei)+^ 2 \Mi\\b^ 

S ' £ 3 = 1,2 

and by H 2 £ the closure of C 00 with respect to the norm 
(2.9) 



+e 



£ I|9 s ^IIl2(v)+ 5 ll^wllraoo 



2 



3=1,2 



We notice that H$ and coincide with 



#L = ^,)€ 



(2.10) 



(2.11) 





=0,1 ^ 2 L=0,1 


= 0inL 2 ((0,oo)); 




0,1 


Oin L 2 ((-l,l)) 




=0 = ^| s= _ i; M s=0 = i>v 


s=l 


^ 2 ((o,i))} 




§ = 


(■01,-02,^) G I 


4e|||*|U 2 <oo; 










=0,1 ~~ ^ 2 L=0,1 


= in L 2 ((0,oo)) 




=0,1 


0inL 2 ((-l,l)) 


^L= 


=0 = ^L=-l ' 


M s=0 = i>v 


1 1 i n 

l s=l 


^ 2 ((o,i)); 




9s 


1 

s=0 £ 5s 


dijj2 
•=-i' 9s 


s=0 


1 

e 9s 


inL 2 ((0,l))l. 

s=l J 



In formulae (|2.1U|) and (|2.11|) . in the boundary values, the symbols \ u _ 1 and , 1 must be in- 
tended as trace operators, see, e.g., [32]. We recall that by standard Sobolev trace theorems, see, 



e.g., [32J, for any < 6 ^ e ^ 1 and p = 1, 2, \P G i?? implies 



J| s =0'^ls=±l 



e^-V2((o,i)), 



$*l«=oi G # p ~ 1/2 ((0,oo)) and ^| u=01 G ^^((-l, 1)), with j = 1,2. Moreover § G #f £ 
implies ^| G #V2((0, with j = 1;2 . 



For p = 1,2, and < a, 6 ^ oo, we shall denote by H p ((a,b)) the standard Sobolev space of 
order p associated to the interval (a, 6), see, e.g., f32] . 

Proposition 2.1. The Hamiltonian Hs t£ is unitarily equivalent to the Hamiltonian Hg i£ ■ 
D(H s ,e) C t £ M e defined by 



(2.12) 



D(H 5 



Hi 



(2.13) 



O 2 



i a 2 



ds 2 5 2 du 2 
where Lgi s denotes the linear operator 



'1, 



O 2 



1 5 



|2 1 



<9s 2 <5 2 9n 2 



(2.14) 

with 

(2.15) 



<5/e :r 



98 Is ds 2 



d_J_ 

ds g s /e 



i d 2 

{6/e) 2 du 2 



Ws/ S = 



I 2 is) 



uS/e ; y(s) 



■xS/ejis))' 



"4(1 + uS/e-f(s)) 2 2(1 + u5/e-f(s)) 3 4(1 + u5/e-i{s)) A 



Proof. We prove first that the quadratic form Qg £ in M.g e is unitarily equivalent to the quadratic 
form Q$ £ in M £ defined by 



poo pi 
Q S ,efo*}--= E / / 

.~Tn JO JO 



fc=1.2 ' 



9s 5s 5 2 du du 



dsdu 



+ e / / — I - - h W 5 / e <t>vVv dsdu 



-1 JO ^ V55/e 5s 9s 



+ e 



1 / 50„ <9V>„ 



+ Ws/ e <f>vil>v jdsdu 



with D(Qs, £ ) = -ffj £ (see formula (|2.10p for the definition of Hg E ) and 



(2.16) 



OS 



i a 



is: 



-1/4, 



3/4 a s ^/ e 



+ 1/2 

Stye V 



(2.17) 



1/4 5 



, 1/2 f 9 _i/4 



Let us denote by £ the restriction of the quadratic form Qg £ to C°° , see the definition (I2.4|) . 
Since g§ i e 6 C°°(y) and the bounds (|2.2p hold true, the norm 



1*1 



I*IIL + Q*[*]) 1/a , 



Qs,t 



is equivalent to the iJj £ -norm defined in equation (|2.8p . Then the quadratic form Q^g is the 
closure of £ in the norm || • • For any function \& 6 C 00 , see the definition (|2.4|) . one has 

* = U^^ty G C 00 , where C/^g is the unitary map defined in (|2.7p and C/^ 1 its inverse. 
Let us define 
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From the definitions of Us e and Q$ £ we have 



/'OO rl 
k=l,2 J ° J ° 



dipk 1 d<p k difj k 



ds ds 5 2 du du 



dsdu 



+ 



1 rl 



-1 J 



1 9 I -V41 \ 5 / -1/4 7 \ . 



£"0 



1/2 Q s ^S/e vv, ds \ys/s 
S/e 



1/2 

%£ 9 _i/4 7 x ^ _i/ 4 7 . 



edsdu . 



We notice that 
»i /-l 

-l Jo 
-1 /•! 



1 ^ / -1/47 s ^ _i/4 r ' 

^T7I^(% £ <^(% £ fa) 

e 



£ ^5/ 



edsdu 



-lio e2 ff<5/£ 9s 5s 



edsdu 



+ 



1 rl 



-1 JO 



1 <9 ~i/4\f ®fa 7 



,2 3/4 



<9.s- 



^ + fa 



ds 



+ ^(I%1 /4 J ^ 

£ 9s/ £ V 



edsdu . 



and 



1 1 1/2 

J ^hto^h**^^ 



i ri 



1 JO 



5 2 du du 



edsdu 



+ 



1 1 r- I/ 4 



— — l/) V + 

an an 



1/2 2 



S 2 V an 



edsdu . 



By integration by parts in s one has 



1 rl 



l Jo 
i rl 



1 d -i/4, / 7 . "T^» 



r 2 3/4 9 s Vy ^ 



ds 



ds 



+ 



d 



»<5/ 



edsdu 



-1 JO 



1 



where W$/ e was defined in equation (|2.16p and we used the fact that the boundary terms are 
null because 7 G Cq°((— 1, 1)). By integration by parts in n one has 



1 1.1/4 

1 '"' r % £ d _i /4 



-1 JO 
1 rl 



1/2 2 

(9s/ J ) ( — ^ + fa-g^ j + 72- ( j <^ 



£2 foVitf/e I ^ 



edsdu 



1 



/ ToWs/efafa^dsdu. 
1 Jo d 



where Ws/ £ was defined in equation (|2. 16[) and we used the fact that the boundary terms are 
null because fa\ u=0 1 = 4> v \ u=0 1 = 0. By noticing that Qf e [$, ^] = Qg t£ [&, SSf] one deduces the 
unitary equivalence of the quadratic forms Qs £ and Qs £ - 



Let iif e the Hamiltonian in M e associated to the quadratic form Qs, £ - By integrating by parts 

in s and u, one sees that for any $,f £ D(Qs £ ) 
(2.18) 



fc=l,2 







/ 

1 Jo 


[?( 


1 _ 


dip! 


&L=o 


ds 



i s 2 v>„ 



9s 2 <5 2 9« 2 
5 1 



dsdu 



9s/s 



+ 



ds 2 
1 <9-0i 



s=0 



e 9s 



9s S«5/ £ 



i>v + {W s/£ + (e/5) 2 W 5/£ )^ 



-i 



I 8=0 



9^2 



ds 



1 



s=0 



e ds 



s=l. 



. 1 ^ 
<5 2 <9tt 2 



edsdu 



where we used the fact that g$/ £ (±l,u) = 1 and the fact that £ Hg £ . A straightforward 
calculation gives 

Then by the first representation theorem, see, e.g., |23j . we have that the Hamiltonian Hg E acts 

as H$ t£ defined in equation (j2. 13[) . Moreover the domain of &f e is given by the functions in Hg £ 
such that the boundary terms in equation (|2.18p are zero and 

ll^t^lk <0 °- 

It is easy to convince oneself that since ggu £ C°°(V) and the bounds (|2.2p hold true, one has 

□ 



D ( H 6,e 



' — Hg s which implies Hf £ = H$, £ - 

3. Approximate solution of the resolvent equation 



In this section we give an approximate solution to the resolvent equation (Hs t£ — z)^> = E, for 
z £ C\M and for some suitable choice of 3 £ H e , see Theorem [1] below. In the analysis we are 
forced to renormalize the spectral parameter z to z + n 2 7r 2 /<5 2 , with n integer. In the final part 
of the section we discuss the behavior of the approximate solution as e — > 0. 
We denote by h v the Hamiltonian in L 2 ((— 1, 1)) 

(3.1) D(h v ) := {y £ H 2 ((-l, 1))| y'(±l) = 0} 



(3.2) 



ds 2 



7 2 (s) 



For any function y £ L 2 ((— 1, 1)) we use the notation 



\y\ 



v,s 



\y\\L 2 ((-i,i)) ; l|y|U,s,oo := ||y|U°°((-i,i)) 



and similarly for the scalar product in L 2 ((— 1, 1)). For n = 1,2,3, ... we denote by y n (s) the 
(real) eigenfunctions of h v and by X n the corresponding eigenvalues arranged in increasing order 

7 2 

KVn = -y'n ~ "j2M = KVn \ y'n{± 1 ) = 11 £ N . 



We assume the normalization (y n , ym)v,s — ^n,mi n,m £ N. For any z £ C\M we denote by r v {z) 
the resolvent of h v , 

r v (z) := (h v - z)~ l ; z£C\R. 
For z £ C\R the operator r v (z) is bounded^, in fact ||^(-z)||.^(z / 2((_ 1>1 ))) ^ 1/| Im z\, and r v (z) : 
L 2 ((— 1, 1)) — > D(h v ). The integral kernel of r v (z) can be written as 

( 3 3 ) „ <„. „ J\ — yn(s)y n (s') 



r v (z;s,s') = ■ 



A,, 



^By SS{-) we denote the Banach space of bounded linear operators in some Hilbert space. 



For n large enough one has (n — l/2) 2 7r 2 /4 < A n < (n + l/2) 2 7r 2 /4 and ||y n ||« lS)00 < c where c 
does not depend on n (see, e.g., [36]). Thus the series in (|3.3p converges absolutely and pointwise 
for s, s' G [—1, 1]. As a function of z, the operator r v (z) can be analytically continued to a linear 
bounded operator for z G C\{A n }. 

We give a formula for the kernel r v (z; s, s') which does not involve series (see, e.g., |22t Ch 4.2]). 
For any z G C\M, let ( v (z) and r] v (z) be two generic solutions of the equations 

- + (-t 2 /4 - *)&(*) = o ; -i) = o 

( " 1 -<(z) + (- 7 2 /4-z)^(z) = 0; »£(*;1) = 

and let W v (z; s) be the Wronskian 

W v (z; s) := rj v (z; s)C' v {z; s) - Q v (z; s)rj' v (z; s) . 

We notice that W v (z; s) does not depend on s, see, e.g., [22], and we set W v {z) = W v (z;s). 
Then the integral kernel of r v {z) reads 

( Cv(z; s)r] v (z; s') 



(3.5) r v (z;s,s') 



W v (z) 
Vv{z; s)( v (z; s') 



s > s 



Wv(z) 

We distinguish two cases: 

Case 1. Zero is not an eigenvalue of the Hamiltonian h v defined in (|3.ip - (|3.2p 

Case 2. Zero is an eigenvalue of the Hamiltonian h v defined in (|3.ip - ()3.2p . Then we denote by 

n* the integer corresponding to the eigenvalue zero, i.e., A ra * = and by y* = y n * the 

corresponding eigenfunction, i.e., 

(3-6) Ky* = --^y* - ^y* = 

with y* (±1) = and ||y*|| U]S = 1 (here and below y* denotes the derivative of y*). We 
define the constants 

(3.7) ai :=y*(-l); a 2 :=y*(l). 

By the definition of h v we have that if y n is an eigenfuction of h v then so is y n . Hence 
we can assume that y* is real, which in turns implies a±, a.2 £ R. 

Let us denote by X n( u ) the functions 

(3.8) Xn(u) V2sm(niru) ; n£N. 

Definition 3.1. For any vector S n = (fiXn, f'2Xn, 0) with fi, f 2 G ^ 2 ((0, oo)) and z G C\M, 
with Im y^z > 0, we denote by *$> £ the vector ^ £ = (^le^e^tie), where the functions ipi, e , 
ifj2, £ and ip v>E are defined by 



[(ro(z)h)(s) + qi ^'] Xn (u) 
[(ro(z)f2)(s) + q 2 , E e i ^ s ] Xn (u) 
e[ii >£ r v {e 2 z; s, -1) + ^2, e r v (e 2 z; s, 1)] Xn (u) , 



(3.9) 

(3.10) &, e (*,i0 

(3.11) ^ Vj£ (s,u) 
and where we set 

(3.12) (r (z)/,) (a) := jf ( J /.(s')^' J Im ^ > 0, j = 1,2. 

77ie constants q\ >£ , q2 e , £i e and £2,2 are /ixed oy £/ie relations 

^ ig . fi,e := (Pi + ; Pi := (ro(«)/l)'(0) 

6,e := (P2 + *V^?2,e) 5 P2 := (roCz)/2)'(0) 
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and 



(3.14) 



\Q2,e 



er v (e 2 z; -1, -1) er v (e 2 z; -1, 1)\ /£ 



er„(e 2 z; 1, -1) er„(e 2 z;l,l) / \£ 



Our main result is expressed in the following Theorem [T] and Theorem [2J the proofs of which 
are postponed to Section [5] and Section H] respectively. 

Theorem 1. For any vector S n = (fiXn, f2Xn, 0) with /i,/2 £ L 2 ((0, oo)) let us take ^ £ as it 
was done in Definition \3.1\ Then ^ e £ D(H$ )£ ), moreover for all z £ C\R i/iere exists Sq > 
smc/i i/iai /or a// < e < £o / or all < 6 ^ e the following estimates hold true: 

Case 1. 

5 „~ , 



n 2 7r 2 



«S c- 



•3/2 



Case 2. 



n 2 7r 2 
c5 2 



, 5 5 

< c — + 



,3/2 1 r 5/2 



where c is a constant which does not depend on e, f\, /a and n. 

Remark 3.2. We notice the following estimate for the resolvent of -H^: if S n = (fiXn, f2Xn, 0) £ 
EU then 



S n ,# £ 



S,e 



n 2 7T 2 



-1 



1 



Im z| 



8,e 



n 2 7T 2 



* P - Sr. 



Theorem 2 (Asymptotic behavior of the solution in the edges). Let us take q\ )E , q2, £ , 6l,e and 
^2, e as it was done in equations (|3.13p - (|3.14p . Then: 

Case 1. 



(3.15) 
Case 2. 
(3.16) 




0(e 





[l + 0( 





(3.17) 




a? + a 



where we denoted by Aq the projector 



(3.18) 



Ct\ OL\Oi2 



An 



a 2 + Ctn \ 2 

1 2 \ai«2 



and A = 1 — Aq 



Here and in the following for all a ^ we denote by 0{e a ) a 2x2 matrix such that ||0(e a )||^(c2) 
ce a , for all < e < £0) where c is a positive constant which does not depend on e, p\ and p2- 
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4. Asymptotic behavior of the solution in the edges (proof of Theorem [2]) 

We devote this section to the proof of Theorem [2j We start with the proof of the following 
proposition: 

Proposition 4.1. For any z £ C\R there exists Eq > such that for all < e < Eq the following 
estimates hold true: 



Case 1. 
(4.1) 

Case 2. 
(4.2) 



sup [\r v (e 2 z;s,s')\] ^ c; 

s,s'6[-l,l] 



sup 

s,s'e[-i,i] 



( 2 l ,y*{s)y*{s') 
r v {e z;s,s) + 



e 2 z 



< c. 



where c is a constant which does not depend on e. 

Proof. We notice that in Case 1 the series Y^ n is convergent and we use the formula (|3,3 
for the integral kernel r v (z;s,s'). Then, in Case 1, 

\yn(s)\\y n (s')\ 



\r v (e z;s,s') 



^ K-e 2 z 



I A„. - e 2 z\ 



Then we use the fact that ||y n |Us,oo ^ c where c does not depend on n, see, e.g., [36J, and the 
fact that for e small enough |A n — e 2 z\ > 2|A n |. We obtain 



r v {e 2 z; s, s')\ < max [||y„||* 8i0O ] oiTT ^ 



c. 



To prove the estimate (|4.2|) we use again the formula (|3.3|) ; we write the integral kernel r v (z; s, s,' 



as 



r v (e 2 z;s,s') 



y*(s)y*(s') v-^ y n (s)y n (s') 



e 2 z 



+ £ 



A„ - £ 2 Z 



We recall that we denoted by n* the integer associated to the eigenvalue zero A n * = 0. The 
estimate (|4.2|> follows from the fact that the series Ylnjtn* r^-r is convergent and by the same 
argument used for the analysis of Case 1. □ 

Remark 4.2. We notice that Proposition 14.11 implies that there exists £q such that for all 
< e < e : 

Case 1. 

\\r v (e 2 z; -,±l)||„ iS ^ 21 sup [\r v (e 2 z; s, s')\] ) < c ; 

\s,s'e[-i,i] / 

y*(s)y*(s') 



Case 2. 



, n , y*y*(±l) 

r v {e 2 z;-,±l) + y " 



e 2 z 



<2( sup 

s,s'e[-i,i] 



r„(e 2;s,s') + 



e 2 z 



< c. 



where c is a constant which does not depend on e. 
We are now ready to give the proof of Theorem [2j 

Proof of Theorem [2J We use the equalities (|3.13|) in equation (|3.14|) and obtain 



'er v {e 2 z\ -1,-1) er v (e 2 z; -1, 1)\ /pi + iy/zq^. 
er v (e 2 z; 1, -1) er v (e 2 z;l,l) J \p 2 + i^fzq %ei 
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We denote by A e the matrix 

A, 



r v (e 2 z; -1,-1) r v (s 2 z;-l,l) 



2„. 



r v (e 2 z;l,-l) r v (e 2 z; 1, 1) 

By Proposition 14.11 for e small enough the matrix (1 — ie-y/zA e ) -1 is well defined , and we get 
the formula 

"'/...- . 

(1 - ie^A £ y 1 eA £ 



'Pi 
\P2j 



Then in Case 1, the estimate f)4. 1[) gives 

\r v (e 2 z;(-iy,(-l) k )\^c j,k = 1,2 

which implies A e = 0(1) and, consequently, the first estimate in equation (|3.15p . The second 
estimate in the same equation follows directly from the definition of £i j£ and £2,6, see formula 

Amp . 

In Case 2 we use the formula 

a i a k , V- Vn{{-l) j )yn{{-l) k ) 



r v (e*z;(-iy, (-!)«) 



e 2 z 



An - e 2 z 



j,k = l,2, 



see the estimate (|4.2p in Proposition 14. II Reminding that y*((— l) k ) = for A; = 1,2, it follows 
that 



eA £ = 

ez 



1 ( af a\a.2 



+e 



,Ql02 



/y^ yn(-l)y^(-l) y n (-l)y„(l) 

Z^inj^n* \ n -e 2 z Z^n^n* A„-e 2 2 

y> yn(i)yn(-i) y> y™(i)yn(i) 



a 2 + a| 

ez 



Ao+eAi, 



where the matrix Ai = C(l) (see Proposition ^. ip . We get the following formula for the operator 
(1 - «£ X /iA e )- 1 eA £ 

r 1 „,2 1 -1 

(1 - ie v / ^A e )" 1 eA £ 



■ a l + a 2 a ■ 2 

2 _ A + £ - IS y/zAx 



a\ + a 2 . 



A + e 2 Ai 



Let A be a 2 x 2 matrix and a some complex constant not equal to zero, moreover assume that 
the matrix (a + A) is invertible, then the formula (a + A)^ 1 = 1/a — (a + ^4) _1 ^4a _1 holds true. 
We use the latter formula to obtain: 



(1 - iey/zAe) l eA e = 

Here we used AqAo = and set 
-.2-1-1 



An 



£ 2 A 2 



\fz 



.a 2 + a 2 



.a 2 + a 2 . 



Aq + e — ie 2 \fzA\ 



i^+s 2 A 2 . 
\fz 



Aq + e — ie 2 ^fzA\ 



Ao 1 [e 2 A, 



.a 2 + a 2 



Aq + e — ie 2 \fzA\ 



-1 



■ Q ? + a 2 a • 2 

', _ A + £ - Z£^VzAi 



[ - ie 2 y/zA 1 



.a 2 + a 2 , 



a 2 + a, 2 ' 


-1 - 


r ^ . 













+ a 2 . 



A + £ 2 Ai 



We remark that for e small enough one has ||A2||^(c 2 ) ^ c, where the constant c does not depend 
on £. By using again the formula (a + A) -1 = 1/a — (a + A) _1 ^4a _1 we get 



— £ 



• a l + a 2 a -2 

-An + £ - IE y[zA\ 



An 



Ar 



a\ + a 2 . 



+ £ 



• a l + a 2 a -2 

1 _ A + £ - %e\fzA\ 



-1 



£ — IE y/zK\ 



An 
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Hence we have 



where A3 is given by 



. , — . _ -| Ao Ao 9 t 

(1 — ieWzA £ ) ek e = i—= — e— ^ ^ + e A 3 

^fz ol\ + a l 2 



e 2 A 3 = e 2 A 2 + e 



Aq + £ — ie 2 A/zAi 



-1 



e — ie 2 y/zK\ 



An 



at + a. 



We notice that for e small enough one has ||A3|| ^( C 2) ^ c, where the constant c does not depend 
on s. Hence we obtain, for e small enough, the following formula for q\ £ and qi E 



(4.3) 



An 



An 



+ £ 2 A 3 



'pi 

\P2; 



which yields the estimate (|3.16f) . The estimate (|3,17p comes from the definition of £i )£ and £2,0 
see equation (j3. 13|) . and from equation ()4.3p . □ 

5. Asymptotic behavior of the solution in the vertex region (proof of 

Theorem [T| 

We devote this section to the proof of Theorem [TJ We first give some preliminary estimates on 
the function f/>„ )£ , see Proposition 15.31 below. 

For any z £ C\R, we denote by ri (z) the resolvent of the Neumann Laplacian in L 2 ((— 1, 1)). 
The integral kernel of ri (z) can be derived from formulas (j3.4j) - (|3.5p by setting 7 = 0; a 
straightforward calculation gives 

( COs(y/z(s + 1)) COS(s/z(s' - 1)) 



(5.1) 



y/zsm(2y/z) 

COS(y/z(s — 1)) COS(y/z(s' + 1)) 

v / ^sin(2- v /i) 



s < s' 



s > s 



Proposition 5.1. For any z € C\R there exists £n > such that for all < e < £n ^ e following 
estimates hold true: 



(5.2) 



A r (0)f e 2 _ n 



< c 



where c is a constant which does not depend on e. Moreover for any g £ L 2 ((— 1, 1)) 

d 



(5.3) 



^ C\\9\ \V,S ! 



where c is a constant which does not depend on e. 
Proof. By a direct computation 



(5.4) 



ds 



■i 0) (e 2 z-;±l) 



2 






r: 



d cos(V e 2 z(s ± 1)) 



ds v / £ 2 isin(2v / £ 2 i) 



ds 



sin(v / £ 2 ^(- ± 1) 



I v,s 



sin(2v £ 2 z) 



For £ small enough, the norms || sin(v / £ 2 -z(- — s arid || sin(v / £ 2 ^(- + l))||u s can be bounded 

by 



(5.5) || sin(v / £ 2 z(- - 1) 

and similarly 
(5.6) 



\v,s 



1 

-1 



giV e 2 z(s— 1) g— j\/£ 2 2(s— 1) 2 



ds ^ C£ 2 



si^v^o + i))!! 2 ^ C £ 2 , 
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where c is a constant which does not depend on e. Using the latter estimates in equations ([57 
we get 



sin(2Ve 2 z)| 2 

which imply estimates (|5.2[) , Also the estimate (|5,3|) can be obtained through a direct calculation. 
For any g £ L 2 ((-l, 1), ds) 



sin(v^(s + l)) f 1 



Then 



(5.7) 



ds 



sin(2y/z) 
sin(y^(g - 1)) 
sin(2\/i) 

sm(Ve^(- + 1)) f 1 



cos(y/~z(s' — l))g(s)ds 
cos(y/~z(s' + l))g(s')ds' . 



sin(2ve 2 z) 
sin(Ve 2 -z(- — 1)) 



sin(2ve 2 z) J—i 

From the Cauchy-Schwarz inequality, for all s 6 ( — 1,1) 



J cos(V e 2 z{s ± l))g(s')ds' 



-i 

cos(v / e 2 z(s / — l))g{s')ds' 
cos(\/e 2 z(s / + l))g(s')ds' 

^ clloll,,.., , 



which implies 
sin(\/e 2 i(- ± 1)) Z" 1 



cos(ve 2 z(s =F l))g(s )ols 



1 



sin(2v e 2 ^) 

Using again the estimates (|5.5p and (|5.6p we obtain 
sin(\/e 2 i(- ± 1)) f 1 



„ iS |sin(2Ve 2 z) 



cos(Ve 2 z(s' =F l))g(s')ds' 



I sinfv e 2 z(- ± 1) 



U.S 3 1),, 



□ 



sin(2v e 2 z 

and the estimate (|5,3p follows from the last estimate and from equation (|5.7p 

Proposition 5.2. For any z S C\R there exists Eq > smc/i i/iai /or < e < £o ^ e following 
estimates hold true: 

Case 1. 

d 



(5.8) 

Case 2. 
(5.9) 



^(eV,-l) 



— r„(e js; •, -1 + = 

ds e l z 



^(eV.l) 



— r„(e z;-,l H = 



< c. 



where c is a constant which does not depend on e. 

Proof. To prove the estimate (|5.8|) we use the well known resolvent identity 

r v (e 2 z; s, -1) = [r v (e 2 z)5(- + 1)] (s) = [[rf\e 2 z) - r^ (e 2 z)(-^ 2 /4)r v (e 2 z))5(- + 1)] 
=ri°\e 2 z; s, -1) - [rf (e 2 z)(- 7 2 /A)r v {e 2 z; •, -1)] (s) . 
Taking the derivative of the latter expression we obtain 



(5.10) 



^V,-i) 



^\e 2 z-;-l) - ^[r(°)( £ 2 .)(- 7 2 /4K( £ 2 z;,-l) 
<c[l + ||^(e 2 2;-,-l)||„ iS ] , 



15 



where we used the triangle inequality, estimates (|5.2p and (|5.3p . and the fact that 7 is bounded. 
As we have proven that in Case 1, \\r v (e 2 z; —1)\\ V:S sj c, see Remark 14.21 we get the first 
estimate in the equation (|5.8|) . the proof of the second one is similar and we omit it. 
To prove the estimates (|5.9p , we rewrite the equation ()5.10p as 

yV(-i)' 



r v (e 2 z;s,-l) =r^(e 2 z;s,-l) 



(5.11) 



r^(e 2 z)(- 7 2 /A) ( r v (e 2 z; ■ , -1) + 



e 2 z 



+ 



e*z 



Then we use the fact that y* is the solution of the problem (|3.6j) and we get 

d 2 y*y*(-iy 



e 2 z 



But 



ri°He 2 z) — u 



d^_ 

'ds 2 '' 



^( £ 2 ,)(-7 2 /4)^ r ^ j 

, d 2 



r^(e 2 ^; s. 



1 
1 



ds' 2 



y*(s')ds 



rf (e 2 z; 



cfs 2 e 2 z 



ids' 2 



^(^s^V)^' = -e 2 z[r(°)(e 2 z)y*]( S ) - y*(s) , 



where we used the fact that the boundary terms are null because 



-V* 
ds 

and we also used the fact that 

d 2 



s=±l 



0; ^P>(^., 



0; 



s'=±l 



ds' 2 



rf (e 2 z; s, s') = -e 2 zr<°> (e 2 z; s, s') - 5(s - s') . 



We remark that ri (z;s, s') = ri°\z; s' , s), see equation (|5.ip . We can rewrite the equation 
(j5TTT]) as 



r v {e 2 z;s,-l) =4°\e 2 z; s, -1) 



r(°)(e 2 z)(-7 2 /4) r,(e 2 z;.,-l) 



e 2 z 



h[ri \e 2 z)y*](s) + y *(s)) y -^ 



or equivalently 



r v (e z;s,-l) +y*(s) 



*,„^*(-l) 



e 2 z 



4 0) ( £ ^)(-7 2 /4) f r^ 2 ,; , -1) + 



= r(°)( £ 2 z;S) _ 1) 

-[r(°)( e 2 z)yV(-l)]( S )- 
Now we can give an estimate of the derivative of the latter equation 



e 2 z 



(s) 



[r v (e 2 z; -,-!)] + 



4 0) (e 2 z;.,-l) 



ifV(-i) 



e 2 ,? 



+ 



r(°)(A)(- 7 2 /4)(r B (A,,-l) + ^P 



+ 



d 
ds 



rf\e 2 z)y*y*(-l) 



< c, 



where we used the triangle inequality, estimates (|5.2p and (|5.3p . and the fact that 7 is bounded. 
Moreover we also used the result stated in Remark 14.21 and the fact that [|2/*|| vs = 1. The proof 
of the second estimate in equation (|5.9p is identical and we omit it. □ 
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Proposition 5.3. Let tp Vt£ be as in equation (|3.1ip . Then there exist e$ > such that for all 
< e < £o the following estimates hold true: 

Case 1 
(5.12) 

Case 2 
(5.13) 



L 2 (v) < ce(|£i, £ | + |$2 |e | 



d_ 
ds 



^ ce{\Ci, E \ + |&, e 



with 



0_ 



L2(y) 



^ce(|ei, £ | + |6,, 



KMi/*(-i) + 6, e y*(i)]xn(«) 



where c is a constant which does not depend on e, f\, $2 and n. 

Proof. The proposition is a direct consequence of the definition of ip V)£ , see equation (|3.11|) , and 
of Remark 14.21 and Proposition 15.21 □ 



The following remark makes sense only in the Case 2: 

Remark 5.4. From Theorem [21 equation (|3,17p . and since y*(— 1) = a.\ and y*(l) = a 2 we 
have that there exists eo > such that for all < e < £o 

l|£ 1)£ y*(-l)+& lE y*(l)| ^c(\ Pl \ + \p 2 \) 

where c is a constant which does not depend on e, p\ and p 2 . Thus implying 

d 

lL2(V) 



(5.14) 



< c(\pi\ + \P2\) ; 



<9s 



^c(\ Pl \ + \p 2 \) 



V(V) 



The following two propositions make preparations for the proof of Theorem [TJ 

Proposition 5.5 (Regularity properties of ^& £ ). F° r an V vector H n = (fiXn, f2Xn, 0) with 
/i)/2 £ £ 2 ((0,oo)) Ze£ us ta&e ^ £ as if was done in Definition ^. 1[ Then ^ £ G D(Hs,e)- 

Proof. From the definition of r«(z), see equation (|3.5p . it follows that the functions 7\,(e 2 z; s, ±1) 
are well defined and that, for any e > 0, r\,(e 2 z; •, ±1) G H 2 {{— 1, 1))- Then for any < 5 ^ e ^ 1 
and from the definition of "ipi >s , ip2,e and ip v ,e, see equations (|3.9p - (|3. 1 1 1) . it follows that, 

It remains to prove that the functions ?/>i )£ , ^2,e an d satisfy the boundary conditions given in 
the definition of D(H$ t£ ), see equation (|2.12p . The boundary conditions in u = 0, 1 are trivially 
satisfied because of the definition of the function Xn(u), see equation (|3.8p . 
Since the costants gi j£ and q 2 ^ £ are related to the constants £i )£ and t; 2t£ through the relation 
(|3.14p . one has that the functions ^/>i )£ , ip 2>£ and t/) V}E satisfy the conditions: 

(5.15) Vv(-1, u) = qi,sXn(u) = Tpi, e (0, u) 

(5.16) tp Vj£ (l,u) = q 2 , £ Xn(u) = V>2, £ (0, u) . 
From the definition of r v (z), see equation (|3.5p . one has 



(5.17) 
(5.18) 



d_ , d c^(-z;g)^(z;i) C(*;s) 

ds rv[z,s,i) ^ c , (z;1) 

_d d 7fr,(z;g)C„(z;-l) r^s) 

ds r^z,s, ij ^ Wv ( z ._y rj' v (z;-l) 
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where we used W v (z) = W v (z; 1) = W v (z; —1) in equations (|5.17p and (|5.18p respectively. Since 
C v (z; — 1) = and n' v (z; 1) = the following relations hold true 



Hence 
(5.19) 



d_ 

ds 



r v (z; s, 1) 



d 



s=l 



i ; 



s=-l 



S r w (*;.,-1) 



s=-l 



s=l 



0. 



ds 




^l )£ (0,tt) = (r (2:)/l) / (0)Xn(w) + iVFgi,eXn(«) = Cl,eXn(t 



1_S 

e <9s 
1 9 



^„ )£ (-l,u) 



(5.20) ^^2, £ (0, n) = (ro(z)/ 2 ) (O)xnN + iVz<l2,sXn(u) = &,eXn{u) = "^^Vt 1 ' n ) ' 
and we have proven that ^ 6 D(Hs £ ). 



□ 



Proposition 5.6. For any couple of constants £i e and ^2,e fei us ta&e E as we did in equation 
(IHTTD . i/ien 
Case i. 



(5.21) 

Case £. 
(5.22) 



<5/£ 



n 2 7r 2 



(o»< 



-e 2 z 



n 2 7r 2 



e 2 z 



<c-[|Ci, e | + |6, e |]e; 
L 2 (y) e 



(5 (5 
< c-[|ei, e | + |6, £ |]e + c- 



L"(V) 









ll^*IL 2 (v) + 




LHV)- 



where L$/ £ was defined in ()2.14j) and c and c are too constants which do not depend on e, fx, 
fi and n. 

Proof. We first notice that, because of the regularity properties of the function ip VtE , see Propo- 
sition for any < 5 ^ e ^ 1 the function L$/ £ ip Vj£ belongs to L 2 (V). Then we notice that by 
our assumption 7 G Cg°((— 1,1)) and from the definition of g§/ £ and W$/ e , see equations (|2.ip 
and (|2.15p respectively, one has the following estimates 



(5.23) 



96/e 



l°°(v) 



6 

<: c - 

e 



^ + T 



5 

<c- ; 



L°°(y) 



5 1 



5s 



L°°(V) 



5 

^ c- . 

e 



From the triangle inequality it follows that 



L 



8/e 



n 2 7T 2 



e 2 z 



1 d 2 7 2 



9s/ l 



ds 2 



E 2 Z 



i>v,e 



L 2 (V) 



+ 



* - ! 4 
where we used the fact that 



IL2(V) 



a 1 




9 t 




5s g 5/e 


L°°(y) 


OS 





By noticing that 



9u 2 



^„ jE = n 2 Tr 2 ip v ^ . 



1 5 2 7 2 _ 2 



96/t 



ds 2 



£ Z 



9S/e 



4 



and using the estimates in equation (|5.23[) we get 



L 



S/e 



n 2 TT 2 



(S/e) 



e 2 z 





5 




< c- 




L 2 (V) £ 





d r 




+ 


OS 


L*(V). 
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In Case 1 the statement follows from Proposition 15.31 In Case 2 one has to use the trivial 
inequalities 

II^> £ IIl 2 (V) ^ ll^fjE ~~ IL 2 (V) + ll^e I 

d ? , 9 



Il 2 (v) 



<9.s 



9s 



+ 



L 2 (V) 



d_ 

ds 



r £ 



L 2 (V) 



and the statement (|5.22|) follows from the estimates ()5.13p and from Remark [5? 
We are now ready to give the proof of Theorem [TJ 



□ 



Proof of Theorem [TJ The fact that $f e € D(H$ £ ) was proved in Proposition 15.51 From the 
definition of the Hamiltonian Hg £ and of the vector ifr e we have that 



n 2 7r 2 



1 



flXn, hXm ~o 



n 2 7r 2 



e 2 z 



where we used 

Then 
(5.24) 



d 2 1 d 2 n 2 7r 2 



ds 2 5 2 du 2 5 2 



S/£ (5/s) 2 

4>j,e = fjXn J = 1,2. 



H, 



8,e 



n 2 Tr 2 



1 



r3/2 



L 



S/e 



n 2 ir 2 



(S/e) 



E 2 Z 



4>v,e 



L 2 (V) 



We notice moreover that from the definition of p\ and p 2 , see equation ()3.13p , one has 
(5-25) \Pj\ <: 4fjh 2 ((o,oo)) ; j = 1,2. 

Using the last estimate in equations (|3.17p and (|5.14j) it follows that 



I^*IIl 2 (v) + 









L 2 (V). 



< c 



1/2 



where we used ||H n [|jj = [ll./illi 2 ((o oo)) H^ 2 Hi 2 ((o do)) J an< ^ WFiere c is a constant which does 
not depend on e, f\, f 2 and n. Then Theorem [T] follows from the last estimate, from equation 
(|5.24p and from Proposition 15.61 □ 



6. Limit operator on the graph 
We denote by Mg the complex Hilbert space 

Mg := L 2 ((0,oo)) 0L 2 ((O,oo)) 
with standard scalar product and norm. In Mg we define the following selfadjoint operators: 
Definition 6.1 (Limit operators on the graph). 

D(h dec ) :={(x u x 2 ) €Bg\x 1 ,x 2 G if 2 ((0, oo)); Xl (0) = x 2 (0) = 0} 
h dec (x 1 ,x 2 ) = (-x'l,-x' 2 r ). 

V x (0) 



(6.1) D(h aia2 ) := Ux 1 ,x 2 ) e Mg Xl ,x 2 G ^((0,oo)); 

(6.2) h a ^(x u x 2 ) = (-x';,-x^), 
where Aq is the projector defined in equation ([3.180 . 



xi (0) 
x 2 (0) 



0; A 



x' 2 (0) 







1!) 



We notice that the operators h dec and h" 1 " 2 belong to the family of operators — A^' described in 
the introduction. In particular the operator h dec is coincides with the Dirichlet (or decoupling) 
Laplacian on the graph and the operator h" 1 " 2 coincides with the weighted Laplacian (for 
N = 2). 

For any z G C\M we denote by r dec (z) and r" 1 " 2 (z) the resolvents of h dec and h" 1 " 2 respectively, 
r dec (z) = (h^-z)' 1 and r aiQ2 (z) = (h" 1 " 2 -z)' 1 . We notice that for any vector (/i,/ 2 ) G Hg 



r dec (z)(/i,/ 2 ) = (roW/i,ro(z)/ 2 ), 
where the operator tq(z) was defined in equation (|3.12p . Moreover 

r Qia2 (^)(/i, / 2 ) = (r (z)/i + gi^, r (z)/ 2 + 92 e^') , 

with 

/Ol \ -a /Vl \ 

Pi = (roWi)'(0);P2 := (r (*)/ 2 )'(0); 

where pi and p 2 are defined accordingly to the definition used in equations (|3.13p . 
Let us define the operator P n : M £ — > Mg 

^n* := ((Xn,V'l)i2 ((0il)) ,(Xn,^2)L2 ((0il)) ) G Hg , 

for any vector § = (i>x,i>2,i'v) G H £ . We denote by P* the adjoint of P n , i.e., P* : Hg ->• H e 
and for any vector (31,32) G Hg 

^(31,32) = (3iXn,32Xn,0) G ir e . 

For any vector (31,32) G Hg and any vector \& = (^l, "02, V^) G H e we have: 

((3l,3 2 ),P n *) He = ((3iXn,32Xn,0),^) ie = (P*(3i,3 2 ),^) 6e . 

Theorem 3. For any vector S n = {fiXni f2Xn> 0); wii/i ||H n ||j| = 1, and for all z G C\R i/iere 
exists Eq > suc/z i/iai /or a// < e < £0 anc ^ / or a// < <5 ^ e the following estimates hold true: 
Case 1. 

A i£ (z-A 2 /i 2 )-r fa WP b 



Case ,2. 



r3/2 1 p 5/2 



where c is a constant which does not depend on e, f\, / 2 and n, and 

4 -1 

iW* - nV/^ 1 ) := ( Hz 



Proof. We give the proof for the Case 2 only. The proof in the Case 1 is identical and we omit 
it. Since the operators R(z — n 2 7r 2 /5 2 ) and r aia2 (z) are bounded by | Im z^ 1 , it is enough to 
prove that for any vector (31, 3 2 ) G Hg one has 



(31,32), r n R s , £ (z - n 2 7T 2 /6 2 ) - r a ^{z 



, 5 5 

^ c — + 



r3/2 



5/2 



The following inequality holds true 
(6.3) / 

1 '(31,52), P„i?5, £ (z-nV/<5 2 ) 



,ai«2 



(*)Pn 



5 

((3i,32),P„fe, £ (^-^V/5 2 )S n -l'J) + ((31,32), 



(z)PjjCj-, 
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where the vector \f e is the one given in Definition 13.11 From Remark 13.21 and from the Theorem 
[I] it follows that in Case 2 

( gi ,g 2 ),F n \R 5 , £ (z-n 2 ir 2 /5 2 )Z 



Bg 



,(91,92), RsA z ~ nV/<5 2 )E n - $ 



< c 



+ 



^3/2 1 £ 5/2 

The second term on the right hand side of equation (j6.3j) can be explicitly written as 



{91,92) 



(91,52), ([gi,6 - ffije*^*', [g 2 , e - g2]e 



< c. 



where we used the result of Theorem [2j equation (|3.16p , and the estimate (|5.25p (we recall that 
we are assuming ||H n || fi£ = 1, then \\fj\\L\(o,oo)) ^ 1, j = 1,2). □ 

7. Conclusions 

We used the unitary map Ug e , see equation (|2.7|) and studied the problem in the Hilbert space 

H e . In this section we discuss our results in the more natural Hilbert space H^. 

In our model the width of the waveguide 5 must be intended as a function of e such that 

< 5(e) < e a , for some big enough, positive constant a. In what follows we shall not mark 

explicitly the fact that 5 depends on e. 

Let use define the functions xg n 



X5,n( U ) := 5 1,2 Xn(u) 



— sin niru 




n 



were the functions Xn were defined in equation ()3.8j) . 

The main idea of our approach is to find an approximate solution for the resolvent equation, 
i.e., for any vector in M. s>e of the form E s>n = (fiX&,n, hxs,n, 0), where fi,f 2 G L 2 ((0,oo)) we 
look for an approximate solution of the equation 



(7.1) 



n 2 TT 2 



z G 



In section [3] we showed an explicit approximate solution which is in the domain of the Hamilton- 
ian H$ £ . The following proposition is a consequence of the fact that U$ £ is unitary and maps 
D(Hg £ ) to D(H$ £ ), the proof follows directly from Theorem [TJ 

Proposition 7.1. For any vector Eg^ n = (fiXs,n, f2Xs,n, 0) with /ij/2 G L 2 ((0,oo)) let us take 
\& e = U^^ £ , where the vector ^f e was given in Definition I3.il and U^ £ is the inverse of the 
unitary map U$ t£ defined in equation (|2.7p . Moreover in Case 1 assume 5(e) < e 3 / 2 and in Case 
2 assume 5(e) < e 5 / 2 . Then * e G D(Hs, e ), and for all z G 



(7.2) 
with 



H, 



<5,£ 



n 2 7r 2 



n + $e 



1$, 



< C £ 



where c £ does not depend on f\, f 2 and n, and c £ 



as e 



0. 



To prove the convergence to an operator on the graph (in norm resolvent sense) we notice that in 
the edges the approximate solution Vl/ £ is factorized in the coordinates s and u. The transverse 
component is of the form X5,n an d the longitudinal component can be written as it is done in 
the following equation (see also Definition 13. ip : 



(7.3) 



4>i,e(s,u) = [(r (z)f 1 )(s) + qi, £ e l ^ s ]xs,n( u ) = Xl,e( S )XS,n( U ) 



■4>2,e( S i U ) = [( r 0(^)/2)(s) 



Q2, £ e 



X8,n( u ) ~ x 2,e( s )XS,n( u ) > 
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We interpret the s dependent parts in ipx >£ and ip2,e, i-e., the functions xi <£ and x 2 , £ , as functions 
on the edges of the graph and consider the vector (xi i£ , X2,e) £ Elg- By definition for all e > 0, 
we have xi j£ , X2,e G ff 2 ((0,oo)) and these functions satisfy the equation 

As we are able to compute the limit of x £< i(0) = qx t£ , £ e ,2(0) = q 2 ,£-, x' £ 1 (0) = £i j£ and x' e 2 (0) = 
^2, e as e — > 0, we can prove that, see Theorem 
in Case 1 

lim x £ i(0) = , lim x £ 2(0) = , lim x' £ -.(0) = p\ , lim x' 2 (0) = p 2 , 

implying that the limit operator on the graph is the Laplacian with Dirichlet conditions in the 
vertex; 

while in Case 2 

(7.4) lim \a 2 x £ i(0) - «ix e o(0)l = , lim \a lX ' £ ^0) + a 2 x' £ 2 (0)1 = 

e— >0 e— s>0 ' 

where ol\ and 02 are two real constants defined by the zero energy eigenvector y* of the Hamil- 
tonian h v , see equations (|3.6p and (j3.7|) . In this case the limit operator on the graph is the 
Laplacian with a weighted Kirchhoff condition in the vertex. In Case 2 one of the constants a\ 
or a 2 may be equal to zero but they cannot be both equal to zero; in this case the limit operator 
on the graph is defined by a Dirichlet condition on one of the edges and a Neumann condition 
on the other. 

We remark that the analysis performed here can be used also to prove the results of [6] and [7] 
which cover some generalizations of the model presented in pp. 

We also notice that our method applies to settings with one or both edges having finite length. 
The form of the solution in the vertex region (i.e. the function ift Vj£ in Definition 13 . 1 1) does not 
change, while the solution in the edge regions must be adapted to fulfill the boundary condition 
in each endpoint. 
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